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A magnetoconductivity formula is presented for the surface states of a magnetically doped topo- 
logical insulator. It reveals a competing effect of weak localization and weak antilocalization in 
quantum transport when an energy gap is opened at the Dirac point by magnetic doping. It is 
found that, while random magnetic scattering always drives the system from the symplectic to the 
unitary class, the gap could induce a crossover from weak antilocalization to weak localization, tun- 
able by the Fermi energy or the gap. This crossover presents a unique feature characterizing the 
surface states of a topological insulator with the gap opened at the Dirac point in the quantum 
diffusion regime. 
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Topological surface states, composed of an odd num- 
ber of massless Dirac cones, are peculiar to three- 
dimensional (3D) topological insulators [H-Q- Electrons 
in these states have a helical spin structure in momen- 
tum space, and acquire a tt Berry's phase after complet- 
ing a closed trajectory adiabatically around the Fermi 
surface. The tt Berry phase could lead to the absence of 
backscatteringQ, weak antilocalization^], and the ab- 
sence of Anderson localization 'i', I?!. In the quantum dif- 
fusion regime (mean free path ^ system size ^ phase 
coherent length), an electron maintains its phase coher- 
ence after being scattered by static centers for many 
times. As a result, the destructive interference due to 
the TT Berry phase can give a quantum enhancement to 
the classical electronic conductivity, leading to weak an- 
tilocahzation (WAL)[i, Applying a magnetic field 
tends to break the destructive interference, giving rise 
to negative magnetoconductivity (MC), a key signature 
of WAL. WAL is expected in systems with symplectic 
symmetry. Much effort has been devoted to observing 
WAL in graphene 

However, graphene has two 
valleys of gapless Dirac cones with opposite chiralities, 
and the intervalley scattering will inevitably suppress 

WAL SEMI. 

In contrast, the surface states of recently 
discovered topological insulators ^2Te3 and Bi2Se3 have 
only one helical Dirac cone 
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and WAL is intrin- 
Many observations of WAL in Bi2Te3 and 
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sic to them 

Bi2Se3 have been reported recently 
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In particu- 



lar, there is great interest in the effect of magnetic dop- 
ing, which is considered to be an efficient way to open 
an energy gap in the Dirac cone by breaking time rever- 
sal symmetry (TRS) j23l - [25j . This gap is expected to 
give rise to many interesting phenomena, such as Majo- 
rana fermion[2^, topological magnetoelectric effect ! 
and quantized anomalous Hall effect 28|. These devel- 



opments call for a thorough theoretical investigation on 
WAL in topological insulators, in particular, in the pres- 



ence of magnetic doping. 
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FIG. 1: (a) A topological insulator with magnetic doping on 
the top surface, (b) Magnetic dopin g m ay open a gap (A) at 
the Dirac point of the surface states [241. |25|. Ep is the Fermi 
energy measured from the Dirac point. A tilts in-plane spin 
polarization of the massless Dirac fermion out of plane, lead- 
ing to the deviation from the tt Berry phase, (c) Scattering 
of an electron by random magnetic impurities. Dashed lines 
represent the trajectory of the electron. 

In this Letter, a MC formula is presented for the mag- 
netically doped surface states of a topological insulator 
[Fig. [I{a)]. We assume that the mean field produced by 
magnetic doping may open a uniform gap at the Dirac 
point [23, Hsj [Fig. [lib)], and the local fluctuation over 
the mean field can scatter conducting electrons in a ran- 
dom fashion [Fig. [Tf c)]- With the help of the diagram- 
matic technique [iji-Il^, [li-lsi] , we obtain the MC for- 
mula, which consists of two competing terms. Besides 
the WAL term due to the gapless Dirac fermion, an ex- 
tra weak localization (WL) term arises as a result of the 
gap opening. We find that either the gap or the magnetic 
scattering can drive MC of the system from WAL to a 
parabolic dependence on the magnetic field (^ B^). Fur- 
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ther increasing the gap/Fermi energy ratio may drive the 
system to WL after reaching the regime. A crossover 
from WAL to WL is thus expected, tunable by the size 
of the TRS-breaking gap or the position of Fermi energy. 
Beyond the theories for non-Dirac systems with strong 
spin-orbit couphng 30] or strong ferromagnetismjs^, the 
observation of the crossover will provide an important 
signature in transport experiments for the existence of a 
TRS-breaking gap in the topological surface states. 

We describe the magnetically doped surface states of a 
3D topological insulator by the massive Dirac model H 
in a random impurity potential U{r). The Hamiltonian 
of the massive Dirac model is given by 



H 



(Jykx) 



are the Pauli matrices. 



(1) 



where CFx,y,z are the Pauli matrices. 7 ~ hvp, with vp 
the Fermi velocity. The ct^ term represents a gap opened 
by breaking TRS, with — ^q^bB + ^, where the first 
term is the Zeeman energy of the out-of-plane magnetic 
field B, with g the g-f actor and the Bohr magne- 
ton. A is the gap opened at the Dirac point 2J, |25| , it 
originates from the mean field produced by the magnetic 
doping. The Hamiltonian describes two energy bands as 
shown in Fig. mb). In this work, we assume that the 
Fermi energy Ep is tuned into the gap of the 3D bulk 
bands, and intersects with the upper band of the surface 
states, which is the only band included in the calculation. 
Its band dispersion is given by ek — \J -I- h?^ and 
the wave function V'k;(r) = [a, — ie*'^6]"^e*'' "'/-\/5, where 
tan (/? = ky/kx, k is the wave vector, a = cos |, 6 = sin |, 
and cos 9 = h^/ \Jh1 + j'^k'^, S is the area. In this work, 
all the physical quantities will be evaluated at the Fermi 
energy Ep at low temperatures. The density of states 
at Ep is Np = Ep The scattering by nonmag- 
netic and magnetic impurities is modeled by the random 
potential 



C/(r) 
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where a runs over 0, x,y, z. uo is the 2x2 unit matrix. 
are the positions of the randomly distributed impurities. 
Mq depicts the potential at R^ for nonmagnetic impurity, 
and w^^y^z for magnetic impurity. Note that ul.^y^^ do not 
represent the total local exchange field produced by the 
impurity at R^, but the local fluctuation over the mean 
field that gives the gap. Therefore, {U{r))i^p — and we 
can still assume the random potential is delta-correlated 
(t/(r)[/(r'))ii„p ^ (5(r — r'), where (...)imp means average 
over impurity configurations, and we follow the practical 
assumption that different types of impurity scattering are 
uncorrelated (lo| . 

The quantum interference correction to conductivity 
of Dirac fermions can be calculated by the diagrams in 
Fig. [21 which are different from those for the usual two- 
dimensional electron gas (2DEG) [ii,[3l|32| in several 
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FIG. 2: The diagrams for the quantum interference correc- 
tion to conductivity of Dirac fermions. The arrowed solid 
and dashed lines represent the Green's functions and impu- 
rity scattering, respectively, (a) The bare [3l|, [s^l and (b) 
two dressed ,10] Hikami boxes give the quantum conductiv- 
ity correction from the maximally crossed diagrams, (c) The 
retarded Green function with the first-order Born approxi- 
mation to the impurity-averaged self-energy, (d) Ladder dia- 
gram vertex correction to velocity 29] . (e) The Bethe-Salpeter 
equation for the vertex of maximally crossed diagrams. 



aspects [13, [2^. (i) Besides the conventional maxi- 
mally crossed diagram (bare Hikami box) in Fig. HJa), 
two dressed Hikami boxes in Fig. [UJb) are also needed; 
each gives —1/4 as the bare Hikami box for the gap- 
less Dirac cone, (ii) The ladder diagram correction to 
the bare velocity = {l/h)de\^/ dk^ [Fig. [2Jd)] must 
be taken into account, which corrects the velocity to 
= 2w^ for the gapless Dirac cone. We generalize these 
conclusions for the gapless case to the gapped Dirac cone 
as follows. 

The arrowed lines in Fig. [2] stand for the impurity- 
averaged retarded (R) and advanced (A) Green's func- 
tions G^^^(w) = l/(w—ek±ifi./2T), where under the first- 
order Born approximation, the impurity-induced self- 
energy is given by the total scattering time r, with 
h/r = 27rX;k'(|t^k'kP)impi5(a; - ek'), where C/k.k' = 
('0k(r)|C/(r)|'0k'(r)) is the scattering amplitude between 
two momenta. It can be separated into \/t — I/t^ + 
I/tiji, where the nonmagnetic elastic scattering time Tg 
is given by h/Te — 2TrNpnoUQ{a'^ + 6^) and the to- 
tal magnetic scattering time can be separated into 
I/tot — 2/tx + 1/tz, with H/tz = 2TTNpnmu1{a^ +b'^) and 
H/tx — 2TrNpnmu'^{2a?b'^). In-plane isotropy {ux = Uy) 
is assumed, uq depicts the average scattering strength for 
nonmagnetic impurities, while Ux.y.z for magnetic impu- 
rities, no and concentrations of nonmagnetic and 
magnetic impurities, respectively. Tg and Tm are related 
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to the elastic scattering length £e and magnetic scattering 
length Ijn by £e = V^Te and £rn = V DTm, respectively. 
D = VpT'^/2 is the diffusion constant. Considering the 



poor surface mobility [18|-|20j, we assume that £e is much 
shorter than the phase coherence length £(/,, as required 
by the quantum diffusion transport. 

In the conductivity diagrams Figs. [2ja) and (b), the 
vertex function F from the maximally crossed diagrams 
usually is proportional to l/g^, where q is the summation 
of momenta before and after scattering. Since F diverges 
as g — > 0, it contributes mainly to backscattering. This 
allows us to sum k and ki for small q first, and write the 
zero-temperature conductivity correction from the bare 
and two dressed Hikami boxes as 



where Q±=Qx±iQy, Q'^=Ql+Ql, Q = vpT sin 9 (q^^qy), 
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and the "Cooperon gaps" go = 2[-2-J 
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1], 



l/r 
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l/r 



- 1]. We note that it is crucial 
to include all the off-diagonal terms of 7 in the calcu- 
lation. Without them, the vertex will be 2 times larger 
(33I 34 1 when going back to the gapless limit [sl], and 
the derived MC formula can not recover to that for 1/4 
of graphene F(q) in Eq. ([3]) can be obtained by 

letting kct = k and k^ = q — k in Fk^k^j , and for 9 — > 0, 
fk. — ^q-k ~ TT. Finally, we collect the most divergent 
terms of the vertex 
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where rjv comes from the correction to velocity from the 
ladder diagrams in Fig. EJd), with 



r(q) 



h/iirNFT) 
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V ' go 52 ' ^ 



51 



(8) 



1 . T 



64 
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and each of the dressed Hikami boxes gives an extra r]H 
contribution as the bare Hikami box, with 



1 / -1 ^ 
VH = -7;{'^-Vv 



(5) 



r]v and jjh reduce to 2 and —1/4, respectively, in the 
absence of magnetic doping [1, [13, [Si . 

In Fig. [2j the total momentum conserves on the in- 
coming and outgoing sides of the vertex F, allowing the 
Bethe-Salpeter equation of the vertex to be written as 
i Tk^k, = F"^,^ + Ek, r°„k,G« G^_k/k,k„ where 
ka -f k^ — q, Icq^^ are the incoming and outgoing mo- 
menta, respectively. For small q, the bare vertex F^^^.^ = 

(t^k;3,k„C/q-k3,q-k„)imp IS fouud aS F^ ^ « 2^1^ + 



-2t:. 



TFT- 



Different from the usual 2DEG, both the bare vertex 
FR t and the advanced Green function G;5 ^ are ex- 
plicit functions of the momentum angle. We propose an 
ansatz to the full vertex function 

h 



F 



kckfl 



2ttNft 



E 



In 



i{nipa-mipfi) 



(6) 



n,me0,l,2 



where 7„m are the expansion coefficients independent of 
(/Jq.,/3- By putting the ansatz into the Bethe-Salpeter 
equation and expanding up to g^, we obtain the 

solution to the expansion coefficients 
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FIG. 3: (a) WL (ao) and WAL (qi) weight factors as func- 
tions of A/2Ef, where A is the gap, Ef is the Fermi energy. 

(b) WL (4) and WAL (^i) lengths as functions of A/2Ef. 

(c) Magnetoconductivity Aa{B) for different A/2Ef in the 
limit of weak magnetic scattering, = 300 nm. £m = 1000 
nm. Ux = Uz is assumed. 



Zero- field conductivity correction cr^(O) can be calcu- 
lated by performing the integral over q in Eq. ([3]) between 
1/^e and 1/^0, respectively In the presence of the 
perpendicular magnetic field B, q^ will be quantized into 
g2 = {n+l/2){AeB/h) = (n-hl/2)/£|, where n labels the 
Landau levels. Summation over n gives the conductivity 
correction a^{B) at finite field 3. The magnetoconduc- 



tivity Acr(S) 
as 



cr^ (B) - cr^ (0) is found for 



> 1 



1=0,1 



^£l 



^2 



(9) 
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with ^' the digamma function, 



Oil 



2(1 + J- + J-)' 
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90 
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and 1/^2 ^ l/£2 4- In the absence of magnetic 

impurities, ao = 0, ai = —1/2, one predicts WAL with 
a prefactor —1/2, consistent with the experimental ob- 
servations 18l420|. For a finite gap, because ao and ai 
have opposite signs, the MC formula has two competing 
contributions, ai leads to WAL, ao to WL. ^o and 
give corrections to i^, in particular, when they are much 
shorter than f^. This formula is the key result of this 
work. 
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FIG. 4: Magnetoconductivity Aa{B) for different magnetic 
scattering lengths £m and A/2Ef- £4, ~ 300 nm. Ua: = Mz is 
assumed. Shorter £m means stronger magnetic scattering. 



We first examine the limit of weak magnetic scatter- 
ing, i.e., e,n > U [l3-i3- We plot MC for different 
A/2Ef in Fig. M^c). For A/2Ef = 0, MC shows a pos- 
itive cusp, which is the signature of WAL. As A/2Ep 
increases, MC gradually develops a B'^ dependence, and 
the system evolves into the unitary regime. Further in- 
creasing A/2Ep will change the sign of MC from negative 
to positive, i.e., a WL-like MC. Different from the usual 
2DEG, the WL-like MC here has a prefactor - 1/2, in- 
stead of 1 31|. This can be seen in Fig. [Sfa), where we 



show the weight factors of the competing WL and WAL 
terms in the MC formula. In the limit of small A/2Ep, 
ai overweighs ao, so MC is mainly contributed by WAL, 
with the maximal prefactor —1/2. In the limit of large 
A/2Ef, ol\ vanishes and ao goes to 1/2, then we have 
WL with the maximal prefactor 1/2. As shown in Fig. 
|3i;b), either h for smaU AI2Ef or ^o for large A/2£'f 
is much larger than l^^ this keeps the system well in- 
side the quantum diffusion regime, and protects WAL or 
WL. For intermediate A/2Ef, where both the WL and 



WAL terms contribute, the weak MC indicates that 
the system is driven from the quantum to classical dif- 
fusion regime due to the effective reduction of by the 
much shorter ^o and ti . The crossover from WAL to WL 
by changing A/2Ef can be understood with the Berry 
phase [1], which is readily evaluated for the surface band 
V'k(r) as 



27T 

I j dip{ V-'k(r) 



^^k(r))^.(l + 4). (11) 



It gives TT for WAL when A/2Ef = 0, and 27r for WL 
when A = 2Ef- A similar argument was also given for 
the gap opened by the finite size effect [37] . In the limit 
of strong magnetic scattering Im <C i<p, both WAL and 
WL are suppressed, as shown in Fig. ID On the other 
hand, because the Fermi energy in the ratio A/Ef can 
be controlled independently by gate voltage [H, [l^ , it is 
possible to observe the transition from negative to posi- 
tive MC by tuning the gate voltage even in this limit. 
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